ARPES and optical renormalizations: phonons or spin fluctuations 
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Improved resolution in both, energy and momentum in ARPES-data has lead to the establishment 
of a definite energy scale in the dressed quasiparticle dispersion relations. The observed structure 
around 80 meV has been taken as evidence for coupling to phonons and has re-focused the debate 
about the mechanism of superconductivity in the cuprates. Here we address the relative merits of 
phonon as opposed to spin fluctuation mechanisms. Both possibilities are consistent with ARPES. 
On the other hand, when the considerations are extended to infrared optical data, a spin fluctuation 
mechanism provides a more natural interpretation of the combined sets of data in Bi2Sr2Ca Cu2 0s+5 
(Bi2212). 



PACS numbers: 74.20.Mn 74.25.Gz 74.72.-h 
I. INTRODUCTION 

With the advent of improved resolution, both in en- 
ergy and momentum, several angular resolved photo 
emission spectroscopy (ARPES) groups have mea- 
sured the renormalized quasiparticle energies and/or 
lifetimeaiiSiMiSiSiLS in the cuprates. The data shows 
that a definite energy scale exists for the renormal- 
ization. This energy scale has tentatively been as- 
signed by someifi' to interactions of the charge carriers 
with phonons. This provocative possibility has again 
brought to the forefront the open question of mecha- 
nism in the cuprates and in particular the possibility that 
the electron-phonon interaction plays an important role. 
While the electron-phonon interaction is widely believed 
to cause superconductivity in conventional materials, the 
gap symmetry in the cuprates is d-waveSiiSiiiiiSii^ rather 
than s-wave which means that it is the projection of the 
electron phonon interaction on the d-channel which en- 
ters the equation for the critical temperature Tc- This 
would imply that some of the details of the electron- 
phonon interaction would be drastically different in the 
oxides compared to conventional metals. While conven- 
tional superconductors do exhibit anisotropjii^ii^ in their 
superconducting gap as a function of momentum k on 
the Fermi surface, the s-channel always overcomes the d- 
channel interaction. Also in dirty s-wave materials with 
the mean free path £ smaller than the coherence length 
^0, the anisotropy is washed out and the gap becomes 
isotropic. On the other hand, in the pure limit there are 
many ways in which this gap anisotropy manifests itself 
in a conventional superconductor lA^ One example is the 
low temperature (T) electronic specific heat which ex- 
hibits the expected exponential dependence exp{— A/T} 
only for temperatures much lower than the lowest gap 
(A) in the system. 

Arguments for why the form of the electron-phonon in- 
teraction can be very different in the oxides as compared 



with conventional metals have been presented in the lit- 
erature. Some of the critical ideas have been reviewed 
by Kulic.^- The main argument for our purpose here is 
that the electron-phonon interaction becomes strongly 
peaked in the forward direction when corrections for 
strong correlations are properly accounted for through 
charge vertices Another related idea is that because 
of the large static dielectric function in the oxides the 
screening is greatly reduced as compared with ordinary 
metals. This also leads to enhanced forward scattering 
as can be seer>i2i2£ from the following simple argument. 
The Fourier transform of the bare Coulomb potential di- 
verges like for small q, here q is the momentum 
transfer. Screening eliminates this singularity and the 
screened potential goes, instead, hke 1 / {q"^ + Xpj,) , where 
Aft is an inverse screening length. For small Aft the 
increased forward scattering has the general tendency of 
increasing the weighting of the higher harmonics in its 
expansion. For example the expansion of a delta func- 
tion (5(q) has equal weight in each spherical harmonics^ 
Strong forward scattering can in fact lead to d-wave su- 
perconductivity which is not suppressed by the strong 
Coulomb repulsion in contrast to the s-wave channeli^iiSS 

While the possibility of phonon induced d-wave super- 
conductivity cannot be eliminated on general grounds, it 
has not been favored in much of the literature. While 
there is no consensus on mechanism, many workers be- 
lieve instead in an electronically driven scenario: for in- 
stance, the t — J model can exhibit superconductivity^ 
A spin fluctuation mechanism, such as is envisaged in 
the Nearly Antiferromagnetic Fermi Liquid (NAFFL) 
model of Pines and coworkersS^iS^ is also electronic in na- 
ture with the exchange of spin fluctuations rather than 
phonons. Another electronic mechanism is the Marginal 
Fermi Liquid (MFL) modcP-^t^L^ which in its original 
form had s-wave symmetry. 

In this paper we attempt to understand the renormal- 
ization effect observed in ARPES within a boson ex- 
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change mechanism. We wih consider exphcitly both, 
phonon and spin fluctuation exchange. We wiU also 
try to reconcile within the same model, both, ARPES 
data on equilibrium quasiparticle properties and the op- 
tical conductivity. Infrared measurements in the cuprates 
have produced a wealth of information on charge dy- 
namics in the oxides^Siffi They have played a key role 
in our present understanding of the microscopic nature 
of the cuprates. In particular, c-axis infrared conductiv- 
ity data have given detailed spectroscopic information on 
the pseudogap^ which has been widely viewed as a di- 
rect manifestation of strong correlation effects. Recently, 
improvements in resolution have also been achieved in 
a6-plane measurements^ 

From reflectance measurements as a function of fre- 
quency it is possible to extract separately the real and 
the imaginary part of the conductivity cr(T, uj) as a func- 
tion of UJ for a fixed temperature T. Within a general- 
ized Drude model, in which the optical scattering rate 
and optical effective mass acquire a frequency depen- 
dence, we can define an isotropic optical scattering time 
Top{T,uj) as T-p^ = (ri^/47r)5ReCT"^(r,w), where is the 
plasma frequency. This quantity can be determined from 
the optical sum rule on the real part of cr(r, w), namely 
duiai{T, u) = (rip/8). The optical scattering rate de- 
fined above is related to the quasiparticle lifetime which 
may be considered to be more fundamental as they di- 
rectly define quasiparticle motion. While these two life- 
times are related, they are by no means identical, how- 
ever, except for a simplified case of isotropic elastic impu- 
rity scattering. Even in this case, if the electronic density 
of states has an important dependence on energy or the 
elastic scattering is anisotropic, quasiparticle and trans- 
port scattering rates are no longer the same. Note that, 
in as much as the in-plane conductivity is isotropic, the 
optical scattering rate is also isotropic and represents a 
weighted average of the more fundamental momentum 
dependent quasiparticle scattering. While there have 
been some attempts in the literature^ to compare both 
sets of data, i.e.: ARPES and optical directly, we will em- 
phasize here that, because they are not so simply related, 
they cannot easily be compared. On the other hand, any 
viable microscopic model needs to be able to provide a 
unifying description of both sets of experiments. In this 
regard, we will emphasize that in addition to the differ- 
ences between quasiparticle and transport scattering rate 
which we have already discussed, they differ in another 
fundamental way. In a boson exchange mechanism the 
charge carrier-exchange boson spectral density (denoted 
by a^F^Lo) in the explicit case of phonons) which de- 
scribes the equilibrium properties can be quite different 
from the corresponding transport spectral density, usu- 
ally denoted by a^^F{u;)m^i^ The origin of this funda- 
mental difference lies in a well-known factor (1 — cos 6*), 
where 9 is the angle between initial and final electron 
momentum undergoing a scattering process. For quasi- 
particle properties such as its lifetime, only the proba- 
bility that an electron of momentum k leaves the state 



|k) is relevant, while for transport, for instance for the 
d.c. resistivity, it is also important where it ends up since 
backward scattering depletes the current more than for- 
ward scattering. This basic physics is accounted for by 
the (1 — cos 9) factor which eliminates forward scattering 
{9 — 0) and weights with a factor two backward scatter- 
ing {9 = it). Serious estimatesiSiiLiS of these two spec- 
tral densities for a model electron-phonon interaction in 
the poor screening regime with explicit inclusion of cor- 
relation effects, have lead to an estimate that af^F{uj) 
may be smaller by about a factor of 1/3 compared with 
a'^F{uj), the corresponding quasiparticle quantity. On 
the other hand, in the NAFFL model the interaction 
with spin fluctuations is believed to strongly peak around 
(tt, tt) in momentum. This means that the scattering po- 
tential itself weights more strongly the backward scat- 
tering processes, and so we expect that the transport 
spectral density will be larger than its quasiparticle coun- 
terpart. The above arguments immediately suggest that 
phonons effects could very well show up prominently in 
quasiparticle properties while at the same time play only 
a minor role in transport, i.e.: in the infrared conductiv- 
ity and vice versa for spin fluctuations. 

In section II we make general remarks about the renor- 
malization effects due to boson exchange mechanisms 
based on perturbation theory. In particular we contrast 
the electron-phonon case with spin fluctuations. In sec- 
tion III we deal with fits to the ARPES data in Bi2212 
and in section IV we consider, in addition, optical con- 
ductivity data. We describe the constraints on micro- 
scopic models that a fit to both sets of data imposes in 
the specific case of Bi2212. In section V we draw conclu- 
sions. 



II. GENERAL RESULTS 

For an electron-phonon system, the renormalization of 
the electronic quasiparticles follows from a knowledge of 
the electron-phonon spectral density a'^F{u!). This func- 
tion which depends only on frequency, and which is lim- 
ited in range to the maximum phonon frequency, contains 
all of the complicated information about electronic band 
wave functions, dispersion relation, phonon dynamics, 
and electron-lattice vibration coupling which is needed 
to compute self energy effects. In lowest order perturba- 
tion theory, the quasiparticle scattering rate T~p^{uj) at 
energy oj is given by^ 

T-p\uj) = 2T: jdua'Fiv). (1) 



For a 5-function Einstein phonon at lje of the form 
a'^F{Ld) = A5{lli - ue) 

('"^ = I ^ < UJE, 
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and we see that loe sets the energy at which T^p{uj) 
jumps from zero to a finite value. For an a^F^oj) dis- 
tributed in energy, the rise will be more gradual but for 
^ > ^Dj with the Debye energy oju, the scattering rate 
will, again, become constant. 

The expression for the optical scattering rate, on the 
other hand, which is also obtained in lowest order per- 
turbation theory is^ 
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dv {u — i')al^F{v) 



(3) 



which gives Tgp{uj) — for w < tj^; in an Einstein model, 
but now for lu > ue 
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which starts at zero with lo — uje and gradually increases 
towards 2t: Atr with increasing lo. This is quite distinct 
from the abrupt jump to 2t:A at u; = loe found for the 
quasiparticle scattering rate. Here the subscript 'tr' de- 
notes transport as opposed to equilibrium properties. It 
is clear that the singularity in t~^{uj) associated with the 
boson energy scale loe is more significant in the quasipar- 
ticle (equilibrium case) than in the transport scattering 
rate. There is yet another difference between quasiparti- 
cle and transport rates. If we denote the electron-phonon 
interaction matrix element by (7k,k',j^ for electrons scat- 
tering from k to k' with a phonon of energy a;,y(k — k') 
(i^ is a branch index) the quasiparticle spectral density 
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F{lj) = J2 ((l5k,k',.|' S{lj - u;,{k - k')) 
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with the double average over the Fermi surface denoted 

by ((•••)) 

J2 Gk,k'^(ek - A*)'5(£k' - 
((G)) = ^ ^ , (6) 

k 

where is the chemical potential, so all integrals are 
pinned at the Fermi surface. The optical spectral weight 
is relatedi^iSi^ instead to 

a^^Fiu) oc ^ ((^\g-kM' (wk - v\^'f S{uj - tj^(k - k')) 

(7) 

The extra factor of Wk — t'k' with Wk the velocity of the 
electron k weights forward scattering by zero and empha- 
sizes most backward collisions. 

The functions a'^F{Lo) of Eq. © and af^F{ui) of Eq. 
are isotropic, momentum averages, given by the right 
hand side of each of these equations. Both involve a dou- 
ble integral over initial and final states |k) and |k'). The 
full complicated momentum dependence of the electron- 
phonon coupling (7k, k',!^ and of the phonon dispersion 



uj^ik. — k') is included with no approximations, although 
after the double summation over momentum, indicated in 
Eq. ©, a'^F{uj) is now isotropic, dependent only on fre- 
quency. This is the function that enters isotropic Eliash- 
berg theory. An anisotropic formulation of the Eliash- 
berg equations^l also exists and in this case a directional 
electron-phonon spectral density distinct for each elec- 
tron state |k) on the Fermi surface, which we denote by 
a^i^(w), would replace its Fermi surface average: 
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To calculate a'^F{uj) from first principles is difficult and 
complexi^ It can, fortunately, be measured directly from 
tunneling data through consideration of current voltage 
characteristics.-22i3S. The resulting function of uj often, but 
not always, looks qualitatively very much like the phonon 
frequency distribution F(w) but this does not imply that 
momentum dependence of coupling and/or of phonons is 
not included in Eq. Q. Any nesting elfects such as are 
present in the nested Fermi liquid theory of Virosztek and 
Ruvalds^S*^ and in the work of Savrasov and Andersen- 
are fully incorporated on the right hand side of Ec^. 
It is only because we take the left hand side from exper- 
iment that these important issues do not become promi- 
nent in our work. 

A further complication arises in the case of coupling 
to spin fluctuations.^"^ In this case the magnetic spin 
susceptibility Xk,k'(t<^) is involved as well as its cou- 
pling to charge. The susceptibility in the cuprates is 
sharply peaked in momentum space at the antiferromag- 
netic wave vector k ^ k' (momentum transfer) for elec- 
tron scattering from |k) to |k'). In this case, what en- 
ters Eliashberg theorji^ is the complex function Xk,k' {^) 
weighted by electron coupling. For instance, in the 
NAFFL model24 



Xk,k'(w) 



XQ 



l+^2(q-Q)2-z- 



q^,qy> 0, (9) 



with q = k — k', XQ the static susceptibility, ^ the mag- 
netic coherence length, Q = (tt, tt) the commensurate an- 
tiferromagnetic momentum, and lommp the characteris- 
tic energy of the spin fluctuations at the point Q. A Fermi 
surface to Fermi surface approximation does not natu- 
rally reveal itself since this may not include transitions 
with k — k' = (tt, tt) for which the susceptibility is the 
largest. However, one can still introduce, as an approx- 
imation, some average effective electron-spin fluctuation 
spectral function denoted by /^x(u;) which can be used in 
the isotropic version of the Eliashberg equations. It rep- 
resents the appropriate average of the susceptibility that 
enters superconductivity. Its value is to be determined 
from consideration of experimental data. As described by 
Carbotte et alm^ this isotropic spectral density is given 
by the second derivative of lu times the optical scatter- 
ing rate T~p^{uj)- RecaU, T~p{uj) = (f7p/47r)3fJeo-"i(w), 
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where (j{u!) is the in-plane optical conductivity. This 
conductivity is isotropic in tetragonal systems. For the 
orthorhombic case, cr(w) with E || a and E || 6 can be 
different. Defining 
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the approximate relation I^^xi^) — W{uj) holdspS We 
emphasize that in terms of microscopic theory W{uj) is 
related to an appropriate average of the full momentum 
and energy spin susceptibility Xk,k'(w) (Eq. ©). Re- 
cently, infrared optical data has been used very effec- 
tively to obtain if^xi'^) in the oxides and we will use 
these in our work:4liJSiiLi& In this way we circumvent all 
of the complications that would arise in a first principle 
calculation of the average effective function Ii^xi'^)- We 
will return to this important point later in equation (|13|) . 
We note in passing that W{lo) very often exhibits a pro- 
nounced resonance like structure particularly if t~j^(ljj) 
data at low temperatures in the superconducting state is 
analyzed. In the following we use the expression 'opti- 
cal resonance' in reference to such a structure. In most 
materials analyzed so far, this optical resonance has a 
spin resonance equivalent measured using inelastic neu- 
tron scattering. We will discuss this in more detail in 
Sec.HVI 

In some sense the opposite extreme to the NAFFL 
model of Pines and coworkers^iiS^ in which the inter- 
action peaks at (tt, tt) is the MFL model of Varma et 
g^j26j27j28 ■y^f}jj(j}j lY^Q interaction with the fluctuation 
spectrum is thought of as momentum independent in a 
first approximation. As originally conceived, this model 
leads to an s-wave superconducting gap rather than d- 
wave as is now generally believed to be the case. This 
represents a real limitation for the model. Nevertheless, 
it has been very useful in correlating much data on the 
normal state transport properties in the cuprates. 

In as much as it is only the magnitude and energy 
dependence of the resulting average spectral function 
a'^F{bj) (/^x(a;)) that matter, these may not be very 
different between the MFL and NAFFL models. Both 
have an average spectral density which is reasonably flat 
as a function of energy and which extends to some high 
frequency of order several hundred meV. Of course, if 
anisotropics on the Fermi surface were to be considered, 
the directional a1_F{uj) {l'^x{^)), Eq. ©, would be ex- 
pected to vary strongly in the NAFFL model and not 
in the MFL model, but here, for simplicity, we have as- 
signed to each electron the same average interaction. For 
phonons in conventional metals the mass renormalization 
parameter Ak = 2 J^j dio a^F [uj) / uj can vary by several 
tens of percentage points over the Fermi surface. For 
the NAFFL model based on the susceptibility l|HJl, the 
variations are much larger, and can be of the order of a 
factor of two or three4^ It is clear that ARPES data in a 
particular direction could considerably under or overes- 
timate the average renormalizations. In this sense optics 
is better. 



We should make one further remark. In conventional 
materials for which a^Fiyj) and a1^F{ijj) have been cal- 
culated from band structure and phonon information 
usually obtained by inelastic neutron scattering, they 
have been found to differii^iiSiSSiSi in shape as a func- 
tion of Lo and in size. But these differences have often 
been overlooked in the literature. In the work of Kulic 
and Zeyher^^'^^'^^ the corresponding differences are large. 
These authors considered directly the cuprates and try to 
account seriously for the ionicity and the increased effect 
of correlations. As described in the introduction, they 
find that the electron-phonon interaction in this case is 
strongly peaked in the forward direction and this leads 
to large differences between a^Fiuj) and a1^F{iS). They 
estimate a'^F{uj) to be about a factor of three larger 
than a1^F{ijj). The consequence of this is that phonons 
would show up much more prominently in quasiparti- 
cle properties (equilibrium) than in the corresponding 
transport property. In particular, at large lo, T^p[uj) and 
Tqp^(cli) (Eqs. l(2Jl and respectively) would approach 
each other if a^Fiyj) and a1^F{bj) were identical, but in 
fact T~j^{uj) is expected to be larger than t~J-{uj). The 
opposite would hold for the NAFFL model which empha- 
sizes backward rather than forward scattering because 
the interaction peaks at (tt, tt). In this case spin fluctua- 
tions should show up more prominently in transport than 
in equilibrium properties. Even though these differences 
are hard to quantify, the general qualitative features of 
transport as compared with quasiparticle spectral weight 
just described must be kept in mind when analyzing data. 
In the MFL model we expect smaller differences between 
quasiparticle and transport spectral densities because the 
underlying interaction is assumed to be approximately 
momentum independent. 

So far we have only discussed scattering rates and have 
emphasized the differences and similarities between in- 
frared optical absorption and quasiparticle lifetimes. In 
ARPES the dressed electronic dispersion relation is mea- 
sured and denoted by fi'k- It is related to the bare band 
dispersion £k through the electron self energy. In the 
electron-phonon case with an Einstein phonon spectrum 
and at T = (zero temperature) the self energy reduces 
to a simple form 



^{uj + ?:o+) = 



In 



lue 



UJ 



— iTT6[uJ — LUe) 



(11) 



In Eq. (|ll|l A is the electron mass enhancement param- 
eter defined as m* /m = 1 + A where m* is the renor- 
malized electron mass at the Fermi surface. In terms of 
a'^F{uj), X = 2j^ dva'^F{v)/v = 2A/loe for the simpli- 
fied Einstein case. The imaginary part of Eq. Hll|l just 
gives back t-1(w) = -2SmS(w + iO+) of Eq. ©. The 
real part gives the renormalized energies as solutions of 
equation i?k = £k + 3fieE(£^k) and for ek — > (on the 
Fermi surface) this gives i?k = £k/ (1 + '^): the quasipar- 
ticle mass is simply renormalized to m*. Note that in 
Eq. (|ll|l the self energy has a logarithmic divergence at 
LO = uje aud this leads to a singular structure in _Ek at 
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that energy as has been investigated in detail by Verga 
et al^ to which the reader is referred for details. Here 
we have chosen to emphasize instead the nature of the 
structure in the imaginary part. It is sufficient to re- 
mark that for a finite distribution of phonon energies in 
a^F^Lu) the logarithmic singularity is moderated as it 
also is when temperature (T) is included. In the general 
case, ?R.eT,{uj + iO+) takes the form 



3?eE(u; + iO+) = - j dva^F[v)^e 



^ \ 2 2ttT 



1 V - UJ 
V' - + ^ 



(12) 



where ipiz) is the digamma function. In lead, the 
electron-phonon spectral density a'^F{uj) is well known 
from tunneling spectroscopjii^i^Si^SiSS as well as from di- 
rect calculations. It is shown as an insert in the top 
frame of Fig.^ The phonons extend to 12 meV and show 
a characteristic peak for transverse and for longitudinal 
branches. In the top frame of this figure we show numeri- 
cal results in the normal state for E\f^ on the Fermi surface 
as a function of the bare band electronic energy up to 
30 meV at which point bare (dotted curve) and renor- 
malized energies are pretty well parallel to each other. 
The solid curve applies to the normal state at T = IK. 
We see structure in the renormalized quasiparticle energy 
i?k ranging up to roughly 10 meV. The structure clearly 
corresponds to the structure in the phonon spectrum (in- 
set). Beyond this range multi-phonon processes are op- 
erative and the renormalization effects are less. Further 
results in the normal state are for T = lOK (dashed) 
and T = 100 K (dash-dotted). We see that by the time 
this last temperature is reached, which is of the order of 
the maximum phonon energy, the thermal effects have 
smeared out much of the structure. Note also, that in all 
three cases the curves go through zero at = and the 
slope of each line out of the origin gives the renormal- 
ized effective mass parameter (1 -f A) at the temperature 
T. For the convenience of the reader we included in the 
bottom frame of Fig. ^ results for the superconducting 
state at T = 1 K (black solid) and at T = 6.5 K (black 
dashed) just below the critical temperature of lead which 
is Tc — 7.2 K. We note that in each case the phonon 
structure is somewhat more pronounced than in the cor- 
responding normal state (shown as gray lines). Also at 
£k ^ the renormalized energy £^k now goes to a finite 
value equal to the superconducting gap A(T) at that tem- 
perature. In an ARPES experiment the gap is seen as a 
shift of the leading edge of the electron spectral density 
A(k, w), downward from the chemical potential level. 

Before moving on to other possible models for the 
renormalization we wish to make an important point 
about how some of the ARPES data has been analyzed 
in the literatureiS*^ It is essential to be aware that the 
bare electron band energies are not known independently 
in ARPES experiments. However, electron momentum 
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FIG. 1: Top frame: the renormalized energies -Ek in lead as 
a function of the bare band energy £k; the renormalization is 
due to the electron phonon interaction characterized by the 
electron-phonon spectral density a'^F{io) shown in the inset. 
The bare dispersion is shown as the dotted line. The other 
curves are for the normal state at T = IK (solid), T = lOK 
(dashed), and T = 100 K (dash-dotted). Bottom frame: same 
as the top frame but for the superconducting state at T = 1 K 
(solid) and T = 6.5 K (dashed) just below Tc = 7.2 K. The 
normal state results are shown as gray lines for comparison. 



change {k — kp) can be measured in the particular di- 
rection of interest. The bare energy is then related to 
this momentum change through the Fermi velocity vp- 
The bare value oi vp denoted v^p is then determined 
from _Ek under the assumption that at ~ 250 meV 
renormalization effects have become small and are neg- 
ligible in a first approximation. This assumption can, 
however, result in a significant underestimate of the true 
mass renormalization A involved as is illustrated in our 
Fig. 121 The solid curve gives vs ek for a case where the 
Einstein frequency defining a^F(uj) is taken to be very 
large, uje — leV, compared with the energies of interest 
in this figure. The curve is almost a perfect straight line 
(there is no structure) with slope giving (1-l-A) with A = 1 
by choice. This is to be compared with the dotted line 
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FIG. 2: The renormalized energy iJk as a function of the bare 
band energy Sk for various models for the electron-boson ex- 
change spectral density. The dotted curve is the bare disper- 
sion relation and is for comparison. The solid curve is for an 
Einstein spectrum with uje = 1 eV and a mass renormaliza- 
tion factor A = 1. The dashed line which shows a slight curva- 
ture is for a model in which the charge carriers are coupled to 
spin fluctuations described by an MMP-form, Eq. 11311 . with 
a spin fluctuation energy Usf ~ 250 meV and a mass renor- 
malization A = 1.1. The dash-dotted line is a straight line 
construction through _Ek ~ 250 meV which is used in ARPES 
experiments to determine the "bare" dispersion as it is not 
independently known. 

which gives the bare dispersion and has slope 1. Apply- 
ing the construction that the "ARPES bare dispersion" 
is a straight line going through the renormalized energy 
at i?k = 250 meV would give the same straight line as the 
solid curve and, therefore, we would conclude that A = 0, 
i.e.: there is no renormalization. As a second example we 
use for the spectral density a form appropriate to a spin 
fluctuation model. The form is obtained directly from 
consideration of experimental data on the infrared con- 
ductivity in the normal state of the cuprates. It is given 
by the experimental form W{uj) of Eq. (|10|) which can be 
adequately represented by a Lorentzian 

(X W{CJ) (X ^ , g(^max - io) . (13) 

(We will refer to this form in the following as MMP- 
form.) In Eq. (|13|) LUgj is a characteristic spin fluctuation 
frequency taken for illustrative purposes 250 meV and 
A = 2 /p'^"="=di^/2x(j/)/i^ = 1.1 with Wmax = 400meV. 
With this spectral density we get the dashed line in Fig. [3 
which is to be compared with both, the dotted line (real 
bare dispersion) and the dash-dotted line (ARPES con- 
structed 'bare' dispersion). It is clear that the ARPES 
construction again drastically underestimates the spin 
fluctuation mass renormalization giving a value of A ^ 
0.1 rather than 1.1. Also, the curve does not show 
any identifiable sharp structure, instead it changes only 



rather gradually. In order to get the correct value of A 
one would need data up to energies higher than the max- 
imum energy involved in the fluctuation spectrum that 
is responsible for the renormalization. For spin fluctua- 
tion theories this energy scale is set by the magnitude of 
J of the t — J model and is high of order several hun- 
dred meV.^^ The t — J model is widely believed to pro- 
vide the appropriate Hamiltonian needed to describe the 
strongly correlated charge carriers in the Cu02-planes in 
the cuprates. In this instance ARPES experiments at 
much higher energies than presently sampled would be 
required to get the full value of A involved. 



III. COMPARISON OF THEORY AND 
EXPERIMENT IN BI2212 

In the interpretation of their ARPES data some ex- 
perimentalists have considered the possibility that the 
renormalizations are due to phonons. We consider here 
only the specific case of Bi2212 for which there also ex- 
ists a recent set of high accuracy optical conductivity 
measurements 1^ The phonon frequency distribution in 
this material has been measured by incoherent inelastic 
neutron scattering by Renker et al.^'^ Such experiments 
give directly F(uj) and a first attempt at a model for an 
electron-phonon interaction spectral density a'^F(uj) can 
be constructed by multiplying the experimental F{u}) by 
a constant adjusted to best fit the ARPES data on the 
renormalized quasiparticle energy i?k as shown in Fig.|31 
The data'^ at T = 130 K are shown in the top frame and 
at T = 70K in the bottom frame as solid squares. On 
the top horizontal scale we show momentum differences 
k — kp in while on the bottom horizontal scale we 
have the value corresponding to our bare band energy 
which implies a theoretical bare energy of 256 meV at 
k-kp = 0.089 A-i. The fit in both cases was obtained 
with a value of A = 1.46. We see that the resulting the- 
oretical results (solid curves) fit reasonably well the ex- 
perimental data. Here, no attempt has been made to get 
a best fit by varying the shape of the underlying spectral 
density a^F{io). The reader is referred to the work by 
Verga et alm^ for an alternative approach to the analysis 
of ARPES data for the material La2-:jSr^Cu04 (LSCO). 
These authors do attempt to get the shape of the spec- 
trum from the data itself and we will return to this issue 
later. Two final comments on Fig. |21 The dotted curve 
which is the ARPES construction for the bare dispersion 
gives at k — kp = 0.089 an £k = 253 meV which is 
nearly the same as the theoretical value of 256 meV. This 
is to be expected since phonon renormalization effects are 
small at 250 meV and beyond. Secondly, we wish to point 
out that the data at T = 70 K is in the superconducting 
state. The formalism needed to treat superconductiv- 
ity is more complicated than that for the normal state 
which we have sketched, but the ideas and concepts are 
the same and no details are included here. The reader 
is referred to an extensive existing literature «22*^ Here it 
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FIG. 3: Top frame: the renormalized energy as a func- 
tion of bare band energy using an electron-phonon spectral 
density a'^F{uj) to describe the couphng. The solid squares 
are the data of Johnson et. ali^ and the solid curve is our fit. 
The full phonon spectral density revealed by inelastic neutron 
scatterings was used with the area under a^F{uj) adjusted 
to get a mass renormalization parameter A = 1.46. In the 
theory a,t k — kp ~ 0.089 ek = 256 meV while the value 
obtained from the ARPES construction is nearly the same, 
253 meV. The curve is for T = 130 K. The dash-dotted curve 
is for a 5-function model for a^F{Lij) at 70meV and A = 0.98. 
The bottom frame applies to the lower temperature T = 70 K 
and is in the superconducting state. 



will be sufficient to point out that the value of the critical 
temperature for superconductivity in a c?-wave supercon- 
ductor is determined most directly by a different function 
that which comes into the mass renormalization chan- 
nel. For an s-wave isotropic superconductor this issue 
does not arise. The quasiparticle a'^F{uj) would also de- 
termine superconductivity and one could ask if ARPES 
data on A are indeed consistent with the observed super- 
conductivity. But in the cuprates the superconducting 
gap has d-wave symmetry. In this case it is the pro- 



the d-wave channel which determines the spectral func- 
tion which is most directly involved in superconductivity 
(rather than its projection on the s-channel). ARPES 
does not measure this projection separately and so this 
technique, stricktly speaking, remains silent on the is- 
sue of mechanism for superconductivity even if phonons 
should be the main cause of renormalization of the elec- 
tron dispersion curves. 

Returning to the top frame of Fig. 01 the dash-dotted 
curve which does not agree well with the data is shown 
for comparison. It is based on a model electron-phonon 
spectral density a^F^uj) in which the coupling is entirely 
to a single frequency fixed at 70 me\i^ with a mass en- 
hancement factor A = 0.98 chosen to get a critical tem- 
perature of Tc = 91 K when s- and c?-channel spectral 
densities are taken to be the same. Such a model does 
not agree well with experiment and when it is used, in 
the superconducting state at T = 70 K, the disagreement 
is even worse. The (5-function model corresponds to an 
extreme in which it is assumed that the electron-phonon 
coupling is to a single mode. This is unlikely to be the 
case. A model based on the entire phonon spectrum is 
more realistic although, as we said before, a^Floj) does 
not need to have the same shape as F(uj). Special modes 
could be weighted more heavily in a'^F{uj) as compared 
to F{u;). 



IV. OPTICAL CONDUCTIVITY 

We can impose a second constraint on renormalization 
effects by considering the infrared conductivity. There 
exists much data on optical conductivity in the cuprates 
and these have given us considerable insight into the 
inelastic scattering involved. The conductivity a{uj) at 
T = (to remain simple) in the normal state also fol- 
lows in a first approximation from the knowledge of the 
electron self energyt^ 



An 



dv ■ 



1 



I](i^) — Yj{lo — v) 



(14) 



jection of the quantity |(7k,k'.iy| 5{uj 



(k - k')) onto 



More complicated formulas exist at any finite tempera- 
ture T and/or for the superconducting case and including 
vertex corrections which we will not repeat here for the 
sake of brevity. They can be found in many places in- 
cluding our own work. Since here we will fit the electron- 
boson exchange density to optical data directly, the re- 
sulting form can be thought of as including higher order 
corrections. 

In Eq. (|14|l is the elastic impurity scattering rate 
while the imaginary part of S(i^) deals with the corre- 
sponding inelastic scattering. What is most often mea- 
sured in optical experiments, is the reflectance R{uj) as 
a function of uj. This is shown in the top frame of Fig.0] 
for Bi2212. The solid curve is the data of Tu et al^ 
for room temperature T = 296 K. The dashed curve was 
obtained with the same neutron based electron-phonon 
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FIG. 4: Top frame: the reflectance for Bi2212 in the infrared 
up to a; = 250 meV. The solid curve presents the data of Tu 
et alw— The dashed curve is for the phonon spectrum and a 
A = 1.46. The temperature is T = 295 K. The bottom frame 
shows the derived optical scattering rate T~p{ui) as a function 
of u). The solid curve presents the data of Tu et al^ The 
dashed curve is a phonon based result which does not give 
enough variation with increasing uj. By contrast, the dash- 
dotted curve fits well and is based on a simple spin fluctuation 
spectrum of MMP-form with uj^f = 82meV and A = 1.79. 



interaction spectrum having a A = 1.46 as we found from 
our ARPES analysis. It provides a poor fit as is illus- 
trated, even more sharply, in the bottom frame of Fig. ^ 
which deals with optical scattering rates. It has now 
become standard procedure for experimentalists to ex- 
tract from their reflectance data the real and imaginary 
part of the optical conductivity and to construct from 
this information the optical scattering rate T^^^ioj). The 
solid line in the bottom frame is the data. The two other 
curves are theoretical and the dashed is based on the neu- 
tron a^F(w) with A = 1.46. It is clear that it does not 
agree with the data, particularly as w increases. This 
defect remains whatever value of A is used. The fun- 



damental problem is that the phonon spectrum cuts off 
too soon as a function of frequency to fit the optics and, 
consequently, the calculated w-dependence of T~p{ixi) is 
much too flat. Any boson mechanism that fits the data 
will need to have an energy scale that extends to order 
of a few 100 me V or so. This is consistent with a spin 
NAFFL-2^ or a MFL modelSS-SLSa but is not compati- 
ble with phonons. It is important to note that this limita- 
tion of a phonon mechanism is independent of any differ- 
ence there might be in the size of a1j,F(uj) as compared 
with a^Fiijj). Both functions cut off at the maximum 
phonon energy (wmax = 80meV in case of Bi2212) and 
this cannot give agreement with the optics. On the other 
hand, a spin fluctuation model provides a natural expla- 
nation because it involves a higher boson energy scale 
of the order J of the t — J-model as previously notediS^ 
This is demonstrated in the lower frame of Fig. ^ by 
the dash-dotted curve which is for an MMP-form with 
LOsf — 82meV and A = 1.79. It flts the experiment al- 
most perfectly with no need for any adjustment of any 
kind. (Deviations from Ugf = 82 meV within ±5 meV 
will change the quality of this fit only marginally.) It is 
interesting to note in closing this discussion that in the 
case of electron-phonon interaction the same value for A 
is required to get a reasonable fit of the ARPES data and 
to get agreement at least in the low energy regime with 
the optical data. This is in conflict with the findings by 
Kufioi^ that in the phonon assisted case a1^F{ijj) should 
be significantly smaller than a^F{ijj). 

In the top frame of Fig. El we compare theory with 
experiment at three temperatures, namely T — 100 K, 
T = 200 K, and T = 295 K. The data on the opti- 
cal scattering rate was obtained by Tu et al^ from 
their infrared measurements in Bi2212 and is denoted 
by gray solid lines (top to bottom is decreasing temper- 
ature). The black dotted lines are theory obtained from 
an MMP-form given by Eq. with ujsf = 82meV. 
The T = 295 K curve repeats the fit shown in the bot- 
tom frame of Fig. 0] The same spectral density, how- 
ever, does not fit well the T = 200K and T = 100 K 
data (dotted lines). To get agreement with these data 
sets it is necessary to augment the MMP-form by adding 
coupling to an additional optical resonance at 43 meV 
(this gives the dashed curves). This is the frequency at 
which inelastic neutron scattering experiments have re- 
vealed a strong peak in the magnetic susceptibility at 
momentum (tt, tt) in the two dimensional Brillouin zone 
of the Cu02 plane;^ in the superconducting state. The 
first report of such a resonance in neutron scattering 
was by J. Rossat-Mignot et al. for YBa2Cu306+5>^ 
We have already described our method for obtaining 
/2,x(w)r^'^^*^Siilii2iSia thus it is not necessary to give 
details here except to mention that Tu et al.^^ also get 
this peak using a slightly different method. We present 
final results for our model Jj^^x([j) in the bottom frame of 
Fig.El While no peak at 43 meV is observed at T 295 K 
(solid line) one shows up at lower temperatures, dot- 
ted line for T = 200 K and a somewhat bigger peak at 
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FIG. 5: Top frame: optical scattering rates for Bi2212 at 
T = 295 K, T = 200 K, and T = 100 K. All are in the normal 
state. The gray solid lines are experiment of Tu et al?^ , dot- 
ted are theory based on an MMP-form for the charge carrier- 
exchange boson spectral density, and the dashed have in addi- 
tion coupling to an optical resonance at 43 meV. The bottom 
frame shows the derived spectral density based on our well 
defined procedure to invert optical data. The T = 295 K 
spectrum has no resonance and is a pure MMP-form with 
LOsj = 82meV (solid line). The other two spectra are mod- 
ified through the addition of coupling to a resonance mode 
at 43meV. Also shown is the phonon spectrum (gray solid 
curve) 1^ 



T = 100 K (dashed curve). At present, neutron scatter- 
ing does not reveal a spin resonance at these tempera- 
tures. Returning to the top frame of Fig. |S1 we stress 
that the introduction of a resonance peak in addition to 
the MMP-form is responsible for the observed rapid rise 
around 50 meV in the optical scattering rate at T = 100 K 
and also, although less obvious, at T = 200 K. The model 
spectral density fits the data remarkably well. Such opti- 
cal resonances have been seen before in YBa2Cu3 06.96'*^ 
as well as in Tl compounds^^*^ in the superconducting 
state of optimally doped samples, and have their equiv- 
alent in spin resonances observed by inelastic neutron 



scattering 



52.53.56.57 



Here the high temperature optical 



resonances in Bi2212 are seen for the first time in the 
normal state above Tc, as reported by Tu et al& Their 
exact microscopic origin, however, cannot be deduced 
from optical data alone. It could be that the 43 meV 
neutron spin resonance forms above Tc in this optimally 
doped material. Alternatively one might argue that this 
may be a phonon contribution (see gray solid curve in 
the bottom frame of Fig. |3J| to the total electron-boson 
transport spectral density. If this were the case, however, 
we would not expect the peak of Fig. |31 (bottom frame) 
to show significant temperature dependence in this tem- 
perature regime. Also, the phonon peak in the bottom 
frame of Fig.|5lhas more weight at low lo than is indicated 
in optics. In terms of the area under the spectral density 
lfj.x{u!) we note that the peak itself contributes less than 
about 15% of the total weight (at T = 100 K), the rest 
comes from the MMP-form. While the issue of the origin 
of the peak is important, here we concentrate instead on 
the consequences of its existence for the ARPES data. 
The possibility of coupling to a resonance mode has been 
discussed in the past in connection with such datai^SiSiffi 

In Fig. we compare our theoretical results for the 
renormalized quasiparticle energies with the ARPES ex- 
perimental data (solid squares) . The theory (solid curve) 
is based on the optics derived spectral density (see Fig. |S1 
bottom frame). The top frame is for T = 130 K in the 
normal state and the bottom frame is for T = 70 K in the 
superconducting state. For the top frame we renormal- 
ized downward the lfj.xiuj) to account for the expecta- 
tion that in a spin fluctuation mechanism the quasiparti- 
cle electron-boson interaction spectral density should be 
smaller than its transport counterpart. To get the fit 
seen in the figure we used A = 1.28. A smaller value of 
A = 1.09 was used in the bottom frame. We see that 
we can get an equally good fit to ARPES with spin fluc- 
tuations including the 43 meV resonance as we did with 
phonons. This new fit (see Fig.EJl, however, has the ad- 
vantage that it can equally well explain the optical data 
while this is impossible with phonons. It is also clear 
from the calculations that the transport electron-boson 
interaction spectral density is larger than the quasiparti- 
cle spectral density by a factor of about two as is expected 
in theories of the NAFFL model. It should be acknowl- 
edged, however, that part of the difference in the value 
of A needed to fit ARPES data as compared with optics 
reflects anisotropics in the quasiparticle rcnormalizations 
over the Brillouin zone. 

Also included in the two frames of Fig. are our results 
for a i?k (heavy dash-dotted) when the resonant peak in 
the spectral density is left out of the calculations. We 
see that the pronounced structure dX k — kp ^ 0.04 
seen in the superconducting state is not reproduced in 
this case. If we used instead the MFL model, i.e.: 
/^x(^) — l'^(^{^max — w), the deviation from the data 
would be worse and even clearly seen in the normal state. 
The resonance peak, absent in the MMP-form and in 
the MFL model, is what gives the observed structure. 
Also shown in the flgure is the the bare dispersion re- 
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FIG. 6: Top frame: the renormalized quasiparticle disper- 
sion _Ek as a function of the bare band energy in meV 
(hght dash-dotted). The sohd curve is based on the spin fluc- 
tuation spectral density with resonance peak obtained from 
conductivity data (Fig. |^ scaled down to a mass renormal- 
ization of A = 1.28, otherwise left unmodified. The heavy 
dash-dotted curve is for comparison and has no resonance 
peak. The light dashed straight line is for comparison and is 
the ARPES "bare dispersion" obtained by making the curve 
go through experiment at 253 meV. This corresponds to a 
theoretical £k ~ 327 meV at — fc^ ~ 0.089 A~^, much 
greater than the ARPES 'bare' value. This frame applies 
to T = 130 K and is in the normal state. The bottom frame 
is for T = 70 K in the superconducting state. 



lation (light dash-dotted curve) on which our calcula- 
tions are based. This is quite different from the ARPES 
derived 'bare' dispersion, the light dashed curve, ob- 
tained by making the curve go through the experimental 
point at 253 meV on the assumption that the renormal- 
ization by boson-exchange interaction has ended at this 
energy. Reference to the difference between the theoret- 
ical bare dispersion (hght dash-dotted) and the dressed 
curve (solid hue) shows clearly that this is not the case 
in a spin fluctuation model. 



It is interesting to note the similarities and differences 
between our work and that of Verga et al^ An important 
difference is that these authors invert the ARPES data 
directly from which they obtain a model electron-boson 
spectral density. Instead, we have inverted the optical 
data which serves as a second constraint on the micro- 
scopies not considered by Verga et al. While they invert 
data in LSCO at three different doping levels x, they, 
nevertheless, find a spectrum which has a peak which is, 
however, not as prominent as either the spin resonance 
peak or the phonon peak in the bottom frame of our 
Fig. O Their peak could simply be a peak of the MMP- 
form as in the solid black curve in the bottom frame of 
Fig. |S1 Nothing definite can be concluded. However, 
we point out that they do find tails at higher energies 
(uj > 100 meV), as we have, which rather support spin 
fluctuations and not phonons. In this paper we consider 
in addition to ARPES optical data and from this we con- 
clude more forcefully that high energy tails extending up 
to several hundred meV exist in the boson assisted spec- 
tral functions describing the electron boson interaction 
in Bi2212. 

Whatever the mechanism by which the quasiparticles 
are renormalized the (ll)-component of the 2x2 Nambu 
Green's function in the superconducting state can be ex- 
pressed in the general formS 



Gii(k,c^) 



(15) 



where the renormalized frequency lu{uj) = luZ(lu). The 
inverse quasiparticle lifetime is deflned as twice the imag- 
inary part of the poles of Gn (k, lu) for a given momentum 
label k. For uj = E'k — iPk we obtain^ 



gk^k2(^k) Aki(^k)Ak2(-Ek) 

^kl(-Ek) E]^Zj 



kl(^k 



(16) 



where the indices 1 and 2 refer to the real and imaginary 
part respectively, and the corresponding renormalized en- 
ergy is a solution of 



AL(i^k) 



(17) 



For the normal state {Ak{u!) — 0) these expressions re- 
duce to those previously given in the limit of weak scat- 
tering with the modiflcation that F is renormalized by a 
mass enhancement factor (1 -I- A). 

In Fig. [7| we emphasize once again the difference ob- 
tained between quasiparticle scattering rates and its op- 
tical counterpart. It applies to T = 80 K for Bi2212 in 
the superconducting state. The gray solid curve repre- 
sents the data of Tu et almL for the optical scattering rate 
and the solid line our theoretical flt to this data which 
also determines the transport spectral density I'^^xi'^)- 
The corresponding value of \tr is 2.18. Consideration of 
the ARPES data on quasiparticle renormalizations shows 
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FIG. 7: Comparison between scattering rates as a function 
of frequency uj obtained from optics and from ARPES. The 
gray solid line is the experimental optical scattering rate at 
T = 80 K obtained from the work of Tu et al.^^ while the black 
solid curve is our theoretical result based on the electron- 
spin fluctuation spectral density Ii^xi^) shown in the bottom 
frame of Fig. |^ (but using data for T = 80 K). The other two 
curves are theoretical results for the imaginary part of the 
quasiparticle self energy. 

us that the equilibrium spectral density must be 

considerably smaller. To calculate the imaginary part of 
T{ijj) of Eq. ifT^ we have reduced if^xi^) by a factor 
of two but without changing its shape. This is done to 
get reasonable agreement with ARPES shown in Fig. |S1 
The results are shown as dashed and dotted lines for the 
antinodal and nodal directions respectively. Not only is 
the shape of the resulting scattering rate different from 
the frequency variation obtained in the optical case but 
its magnitude is also different reflecting the difference in 
magnitude of I'^^.xi'-^) and /^x(^)- 

V. CONCLUSION 

We have reconsidered the interpretation of the quasi- 
particle renormalizations observed in high resolution 
ARPES data for the specific case of Bi2212. We con- 
firm that the observed energy scale is compatible with 
coupling to phonons. We find, however, that present 
data and, in particular, the way they are analyzed might 
well significantly underestimate renormalization due to 
other boson exchange mechanisms distinct from phonon 
exchange and for which the energy scale is larger. As an 
example this would be the case in a spin fiuctuation ex- 
change mechanism as is envisaged in the NAFFL model 
where the energy scale for the spin fluctuations is set by 
J of the t — J model and is consequently very high. 

If in addition to ARPES, optical data is also consid- 
ered, it is found that an explanation in terms of phonons 



becomes less tenable. Of course, there are many com- 
plications and these make a definitive interpretation im- 
possible. What is clear, however, is that transport prop- 
erties and in particular the observed infrared conductiv- 
ity cannot be understood as due to the electron-phonon 
interaction. This holds even if modifications arc intro- 
duced to account for the reduced screening present in 
the cuprates because of the low electron density and also 
the increased effects of correlations. The main point is 
that the underlying charge carrier-phonon spectral den- 
sity is cut off at the maximum phonon frequency and this 
feature is in conflict with optical data. By contrast, spin 
fluctuations give a natural explanation of the data. Our 
analysis of the data of Tu et al^ also provides a natural 
explanation of the ARPES data. What we flnd is that an 
optical resonance is present in the infrared data even in 
the normal state above Tc and this resonance can explain 
the energy scale seen in ARPES. (Below Tc this optical 
resonance has a spin resonance equivalent observed by 
Fong et alm^ using inelastic neutron scattering.) The mi- 
croscopic nature of the optical resonance is not entirely 
known. It could have its origin in spin fluctuations the- 
ories and be the spin one resonance observed in inelastic 
polarized neutron scattering experiments, or could even 
contain a phonon contribution over and beyond the spin 
fluctuation contribution. 

Our analysis of both ARPES and optical data together 
has served to emphasize once more the differences be- 
tween the equilibrium electron-boson exchange spectral 
density a^F(uj) and its transport counterpart af^F{uj) 
be it for phonon exchange or spin fluctuations. The data 
indicates that these two quantities differ in magnitude by 
roughly a factor of two or so. This is understood to arise 
from the fact that transport scattering rates emphasize 
more strongly backward collisions than do equilibrium 
rates. While detailed and elaborate calculations beyond 
the scope of this paper would be required to pin down the 
expected differences quantitatively, present ideas about 
reduced screening in the cuprates lead directly to the 
expectation that phonons will be much less important 
in transport than for equilibrium properties while ideas 
about the NAFFL model lead to opposite expectations. 
Spin fluctuations dominate the transport while at the 
same time could be less important in determining quasi- 
particle properties. This agrees well with the available 
combined set of optical and ARPES data in Bi2212. 
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